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In the present follow-up article of a previous one [1] we illustrate 
the use of the Unconventional Statistical Mechanics described and dis- 
cussed in the latter. This is done via the analysis, resorting to Renyi's 
approach, of experimental results in the case of so-called "anoma- 
lous" luminescence in nanometric quantum wells in semiconductor 
heterostructures, and the so-called "anomalous" cyclic voltammetry 
in fractal-like electrodes in microbatteries. Also a purely theoretical 
analysis is done in the cases of an ideal gas and of radiation comparing 
the conventional and unconventional approaches. In all of these sit- 
uations it is discussed which is the failure to satisfy Fisher's Criteria 
of Efficiency and/or Sufficiency thus requiring to resort to the uncon- 
ventional approach, and what determines the value of the infoentropic 
index in each case, and its dependence on the system characteristics. 
Moreover, on the basis of the results we obtain, it is conjectured that 
the infoentropic index may satisfy what we call a law defining a "path 
to sufficiency". 
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1. INTRODUCTION 

In a preceding article [1] (heretofore referred to as I) has been presented the con- 
struction of a so-called Unconventional Statistical Mechanics (USM), that is to 
say auxiliary forms for use instead of the conventional one - the latter based on 
the quite general and well established formalism of Boltzmann and Gibbs - when, 
as there noticed, the researcher is unable to satisfy Fisher's Criteria of Efficiency 
and/or Sufficiency in Statistics, namely, a failure in the characterization of the 
system and its dynamics in what is relevant for the problem in hands. In the 
present paper we apply the theory to some particular cases, namely two experi- 
mental situations involving systems with fractal-like structures, and an analysis 
of ideal gases. 

We consider in next section measurements of "anomalous" luminescence in 
nanometric quantum wells in semiconductor heterostructures. In section 3 we 
describe the case of experiments of cyclic voltammetry in microbatteries with thin- 
film fractal electrodes, where it is involved the so-called "anomalous" diffusion of 
charges. In section 4 is presented a study of ideal gases comparing the case when 
a proper characterization is used and the one when as incomplete characterization 
is used. In all cases we resort to the use of Renyi statistics [2, 3] described in I; as 
expected the infoentropic index a can only be 1 in the first situation (the criteria 
of efficiency and sufficiency are satisfied), but different from 1 in the incomplete 
description when the criteria of efficiemcy and/or sufficiency are faltering. As 
already noticed in I, the infoentropic index (when different from 1) depends on 
several characteristics of the system and its dynamics, viz. the fractal topography, 
the type and size of the system's geometry, its equilibrium or nonequilibrium 
thermodynamic state, the experimental protocol, and so on. Furthermore, as we 
shall see below, it appears that the infoentropic index can be related to quantities 
characterizing the system and its description, through relations which can be 
considered as indicating a kind of "path to sufficiency" . 

After the presentation of the five sections with the applications, in a last section 
we add some additional comments and a summary of the results. 

2. "ANOMALOUS" LUMINESCENCE 

There exists nowadays a large interest on the question of optical properties of 
quantum wells in semiconductor heterostructures, which have been extensively 
investigated in the last decades as they have large relevance for the high perfor- 



mance of electronic and optoelectronic devices (see for example Ref. [4]). To deal 
with these kind of systems, because of the constrained geometry that they present 
(where phenomena develop in nanometer scales) the researcher has to face diffi- 
culties with the theoretical analysis. A most relevant question to be dealt with 
is the one related to the interface roughness, usually having a kind of fractal-like 
structure, that is, it is present a spatially varying confinement, which leads to 
energies and wavefunction depending on boundary conditions which need account 
for spatial correlations. As a consequence the different physical properties of these 
systems appear as, say, "anomalous" when the results are compared with those 
that are observed in bulk materials. A particular case is the one of photolumi- 
nescence which we briefly describe here. The conventional treatment via the well 
established Boltzmann-Gibbs formalism has its application impaired because of 
the spatial correlations resulting from the spatially varying confinement (as no- 
ticed above), relevant in the characterization of the system, on which one does 
not have access to (obviously the interface roughness varies from sample to sample 
and one does not have any easy possibility to determine the topography of the 
interface). This is then the reason why the criterion of sufficiency is not satisfied 
in this case. 

Let us consider a system of carriers (electrons and holes) produced, in the 
quantum well of a heterostructure, by a laser pulse. They are out of equilib- 
rium and their nonequilibrium macroscopic state can be described in terms of an 
informational-based statistical thermodynamics [5]. It is characterized by the time 
evolving quantities energy and density or alternatively by the intensive nonequi- 
librium variables (Lagrange multipliers in the variational approach to statistical 
mechanics) quasitemperature T* (t) and quasi-chemical potentials /i e (t) and \x h (t) 
(e for electrons and h for holes) [6]. Electrons and holes do recombine producing 
a luminescence spectrum which, we recall, is theoretically expressed as 

/ (U | t) OC /nkx W (0 * " 4k ± " & J , (1) 

where f e and f h are the populations of electrons and of holes, Ml = YvjJ — Eg, 
with uj being the frequency of the emitted photon and E G the energy gap, and the 
e nk]_ are the electron (hole) individual energy levels in the quantum well (index n 
for the discrete levels and for the free movement in the x — y plane) . 

The textbook expression for the energy levels corresponding to the use of 



perfectly smooth bidimentional boundaries is given by 
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where Lqvk is the quantum-well width and w*^) is the effective mass and the 
populations / take a form that resembles instantaneous in time Fermi-Dirac dis- 
tributions [5,6]. Using this expression in the calculation of I of Eq. (1) and 
comparing it with the experimental results [7] one finds a disagreement, and then 
it is used the name of "anomalous" luminescence for these experimental results. 
This is a consequence that we are using an improper description of the carriers' 
energy levels - we are not satisfying the criterion of sufficiency (as discussed in 
I) -, resulting of ignoring the roughness of the boundaries (with self-affine fractal 
structure [8]) which needs be taken into account in these nanometric-scale geome- 
tries, and then the boundary conditions to be placed on the wavefunctions are 
space dependent. Hence complicated space correlations are to be introduced, but 
to which we do not have access (information), as already noticed. Hence, this 
limitation on the part of the researcher breaks the sufficiency criterion, and ap- 
plication of the Boltzmann-Gibbs-Shannon-Jaynes construction is impaired. As 
noticed in the preceding article, one can try to circumvent the difficulty (which, we 
stress once again, resides in the limitations the researcher has to possess a proper 
characterization of the system and its dynamics, and not in Boltzmann-Gibbs 
statistics) introducing unconventional statistics based on parameter-dependent 
structural informational entropies. 

To deal with the "anomalous" luminescence in nanometric quantum wells in 
semiconductor heterostructures, we have used the unconventional statistics that 
is derived from Renyi structural entropy, which depends on a single parameter, 
namely, the infoentropic index a (see I). Taking as basic variables, for describing 
the nonequilibrium thermodynamic state of the "hot" carrier system, the energy 
E c (t) and electron and hole particle number (or density), iV e (t) and jVj, (t) [6] (the 
thermodynamically conjugated intensive variables are the quasitemperature T* (t) 
and quasi-chemical potentials /i e (t) and /i e (t)), and using Renyi informational 
entropy, as shown in I [cf. Eq.(I-42)], we do have for the carriers' populations, to 
be used in Eq. (1), that 

fl h L a (t) = {{l + (a - 1) B*» (t) [et[ ~ tfP (*)] }^ ± I}" ■ (3) 

Here E a and fi a are modified forms of the Lagrange multipliers associated to the 
basic variables energy and particle number (see I). The first, B , plays the role of 



a reciprocal of a pseudotemperature, and /i of a quasi-chemical potential (see I 
and [6]). Using these populations in the nondegenerate limit (cf. Eq. (53) in I), 
which depend on e^|, that is, the ideal single-carrier energy level of Eq. (1), the 
luminescence spectrum is given by 
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where & //a (t) = ^ (t) + £B* (t), m" 1 = [m 
The second contribution in the last term in Eq. (4) is much smaller than the first, 
as verified a posteriori, and then we approximately have that 
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The experiments reported in [7] are time integrated, that is, the spectrum is given 
by 
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where At is the resolution time of the spectrometer. Using Eq. (5) in Eq. (6), 
and in the spirit of the mean-value theorem of calculus we write 
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(7) 



introducing the mean value /3 e ff a (as an open parameter), which we rewrite as 

[<^e//J 1 = ksQa, defining an average, over the resolution time At, effective 
temperature of the nonequilibrium carriers, that is, a measure of their average 
kinetic energy (see [9]). 

In Fig.l is shown the fitting of the experimental data with the theoretical 
curve as obtained from Eq. (7). It contains the results referring to four samples 
having different values of the quantum well width. The corresponding values of 
a, and the kinetic temperature G a are given in Table I. The information-entropic 
index a depends, as expected, on the dimensions of the system: as the width of 
the quantum well increases the values of a keep increasing and tending to 1. This 
is a clear consequence that the fractal-like granulation of the boundary surface 



becomes less and less relevant, for influencing the outcome of the phenomenon, as 
the width of the quantum well falls outside the nanometer scale, and is approached 
the situation of a normal bulk sample. On the other hand the kinetic temperature 
of the carriers is smaller with increasing quantum well width, as also expected once 
the relaxation processes, mainly as a result of the interaction with the phonon 
system, become more effective and the cooling down of the hot carriers proceeds 
more rapidly. 

Moreover, we can empirically derive what we term as a law of "path to suffi- 
ciency" , namely, 

»< £ > = rr£ • < 8 > 

where, by best fitting, L x ^ 139 ± 17 L 2 ~ 204 ± 24, all values given in Angstrom. 
We do have here that as L largely increases, the entropic index tends to 1, when 
one recovers the expressions for the populations in the conventional situation (see 
I), but as L decreases a tends to a finite value Li/L 2 , in this case ~ 0.7 ± 0.06. 
This indicates that the insufficiency of description when using Eq. (2) in the 
calculations (the ideal energy levels) becomes less and less relevant as the size of 
the system increases as already commented. 

It is relevant to notice that there exists computer-modelled experiments, in 
which a certain controlled roughness of the quantum-well boundaries is intro- 
duced, Schrodinger equation is solved and the corresponding energy levels are 
obtained [10]. The conventional statistics is applied, and the results, for this 
model, qualitatively agree with the experimental ones in real systems. Further- 
more, the results of the computer modelled system, with the sufficiency criterion 
being satisfied, can be reproduced using the ideal eigenvalues of Eq.(2) but in 
the Renyi statistics we used adjusting the infoentropic index a [7]. Another ob- 
servation is that experiments in which are present quantum wells with the same 
width, but obtained with careful and improved methods of growing, show that 
as smoother and smoother the boundary surfaces the entropic index a increases 
tending to 1, the value corresponding to the ideal situation of perfectly smooth 
surfaces [7]. Moreover, from the experimental data [7] it can be noticed that the 
linewidth T, in each of the four samples of Fig. 1, appears to depend on temper- 
ature T, for T > 20K, through a power law, say, T ~ T", with v approaching 1 
(the conventional result) as the quantum-well width increases, as expected; study 
of the connection of v and the infoentropic index a is under way. 

As a consequence, it can be noticed that through measurements of lumines- 
cence it is possible to obtain an evaluation of the microroughness of the samples 



being grown, what implies in a kind of method for quality control. On this we 
may comment that the interface structural properties in quantum wells (QW's) 
have been extensively investigated, as they are extremely important for the high 
performance of electronic and optoelectronic QW-based devices[ll-13]. Semicon- 
ductor heterostructures interfaces have been investigated by means of direct or in- 
direct characterization techniques. Direct investigations of the interfacial quality 
have been obtained, for instance, by scanning tunneling microscopy, atomic force 
microscopy, and transmission electron microscopy [14]. However, interfaces are 
not easily accessible through these direct investigation methods, so optical tech- 
niques (which indirectly probe the interfaces) can constitute a useful approach in 
semiconductor-interface characterization in heterostructures, as shown above. 

This first illustration of the theory clearly evidences the already stated fact that 
the infoentropic index a is not a universal one for a given system, but it depends 
on the knowledge of the correct dynamics, the geometry and size including the 
characteristics and influence of the boundary conditions (e.g. the fractality in the 
present case) , the macroscopic (thermodynamic) state of the system of equilibrium 
or nonequilibrium conditions, and the experimental protocol. 

The case we presented consisted of experiments in time-integrated optical spec- 
troscopy. The phenomenon of "anomalous" luminescence in nanometric quan- 
tum wells in semiconductor heterostructures, is also present in the case of time- 
resolved experiments (nanosecond time resolution where the infoentropic index 
and the (nonequilibrium) kinetic temperature change in time accompanying the 
irreversible evolution of the system [15]). The use of USM in a completely analo- 
gous way as done above (note that Eq. (4) is valid for a time-resolved situation) 
allows to determine the evolution in time of the kinetic temperature Q (t) , and 
the infoentropic index a(t), which is then changing in time as it proceeds the 
evolution of the irreversible processes in the nonequilibrium thermodynamic state 
of the carriers, which is a work under way. 

3. "ANOMALOUS" DIFFUSION 

As a second illustration we consider the question of behavior of fractal electrodes 
in microbatteries. As a consequence of the nowadays large interest associated to 
the development of microbatteries (see for example Ref. [16]), the study of growth, 
annealing, and surface morphology of thin-film depositions used in cathodes, has 
acquired particular relevance. These kind of systems are characterized by micror- 
oughneessed surface boundaries in a geometrically constrained region (nanometric 



thin films), and then fractal characteristics can be expected to greatly influence 
the physical properties [8]. In other words the dynamics or hydrodynamics in- 
volved in the functioning of such devices can be expected to be governed by some 
type or other of scaling laws . Such characteristics have been experimentally ev- 
idenced and the scaling laws determined. Researchers have resorted to the use 
of several experimental techniques: atomic force microscopy allows to obtain the 
detailed topography of the surface - and then measuring the fractality of it 
cyclic voltammetry is an electrochemical technique used for the study of several 
phenomena, also allowing for the characterization of the fractal characteristics 
(scaling laws) of the system [17-19]. Particularly, through cyclic voltammetry it 
can be put into evidence the property of the so-called "anomalous" diffusion in 
these systems, what we consider here. 

The difference of chemical potential between an anode and a cathode with a 
thin film (nanometric fractal surface) of, say, nickel oxides, produces a movement 
of charges in the electrolyte from the former to the latter. In a cyclic voltammetry 
experiment these charges circulate as a result of the application of an electric field, 
e (t), with particular characteristics: It keeps increasing linearly in time as e Q + vt, 
where v is a scanning velocity, during an interval, say At, and next decreases with 
the same scanning velocity, i.e. e D + vAt — vt, until recovering the value e Q . A 
current i (t) is produced in the closed circuit, which following the field e (t) keeps 
increasing up to a peak value i p , and next decreases. This current is the result 
of the movement of the charges that keep arriving to the thin film fractal-like 
cathode. 

It is found that there follows a power law relation between the peak value, 
i p , of the current and the rate of change, v, of the electric field, namely i p ~ 
v^. Such current in the fractal electrode, generated by the application of the 
external potential, depends on the charges that are brought to the interface. Such 
density of charges n (x, t) at the interface is the one accumulated as a result of a 
diffusive motion, generated by the difference of chemical potentials between both 
electrodes. Using the standard Fick's law it follows that the power index £ must 
be 1/2. Observation shows that it departs from that value and to account for the 
disparity it has been postulated an "anomalous" diffusion law of the type 

^(r,f)-J} 7 VV(M) = , (9) 

where when 7 = 1 is recovered the standard result. Such kind of result can 
be derived using the USM of I, and we resort to the statistics based on Renyi 



structural information-entropy, as done in the previous section. The question now 
is which is the source of the lack of sufficiency in the well established Boltzmann- 
Gibbs formalism, which forces us to resort to the unconventional approaches. The 
answer resides in that is being used a quite incomplete hydrodynamic approach 
to describe the motion of the fluid. Fick's law is an approximation which gives 
good results under stringent conditions imposed on the movement (one being the 
limit of very large wavelengths, that is, the movement is characterized by being 
of smooth variation in space (near uniform), see Appendix A). Motion of the 
particles (charges) towards the thin film in the cathode proceeds through the 
microroughneessed fractal region, and then involves a description requiring to 
consider intermediate to short wavelengths. Hence, for its description one needs 
to remove the limitations that restrict the movement to be purely diffusive, that 
is, to introduce a higher-order generalized hydrodynamics [21, 22] (see also section 
4). This means to introduce as basic hydrodynamic variables not only the density, 
n (r, t), but its fluxes of all order, I n (r, t) and 1$ (r, t), where r > 2 indicates the 
order of the flux and its tensorial rank: see Eqs. (38) to (42) in Section 4. The 
motion is then determined by a complicated set of equations of evolutions of the 
type [5,21,23] 

p£Hr,t) + V-lM(T,t)=jP(T,t) , (10) 

where r = for the density, r = 1 for the first (vectorial) flux, or current, and 
r > 2 for the all higher-order fluxes, Jn^ are the collision operators, and V- is 
the operator indicating to take the divergence of the tensor. Solving this set of 
coupled equations of evolution is a formidable, almost unmanageable, task. As a 
rule one uses, depending on each experimental situation, a truncation in the set of 
equations (i.e. they are considered from r = up to a certain value, say n, of the 
order r). Furthermore, to build and solve the set of equations with the spatially 
quite complicated interface boundary conditions is practically not possible (at 
most can be attempted in a computer-modelled system). Hence, introducing 
such truncation produces a failure of the criterion of sufficiency when using the 
conventional, and universal, approach. Consequently, a way to circumvent the 
problem is, as done in the previous section, to make calculations in terms of 
unconventional statistical mechanics, what we do resorting to Renyi's approach. 

We take a truncated description introducing as basic variables the densities 
of energy and particles and the first (vectorial) particle flux, the corresponding 



dynamical operators being 
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The unconventional statistical operator is in this case given by Eq. (1.21) 
using in it the auxiliary ("instantaneous frozen" or quasiequilibrium) statistical 
operator 
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where F n (r, t) and F (r, t) are modified Lagrange multipliers that the variational 
method introduces, (3 = 1/ksT with T being the temperature of the system, 
and then we are assuming that the material motion does not affect the thermal 
equilibrium, and fj a (t) ensures the normalization of the distribution, T is a point 
in phase space and we shall designate by dT the element of volume in phase space 
(taken adimensional). 

The equations of evolution for the basic variables, derived in the context of 
the kinetic theory based on the Nonequilibrium Statistical Ensemble Formalism 
(NESEF) [1,5,24] are 



JU(r,0 = -V-I n (r,t) 

(r, t) = -V • /I 2 ] (r, t) + 3 na (r, t) , 



(16) 
(17) 



where the second order flux is given in Eq. (A. 18) of Appendix A and V- is the 
tensorial divergence operator. Equation (16) is the conservation equation for the 
density; the terms with the presence of the divergence operator arise out of the 
contribution resulting from performing the, in this classical case, Poisson bracket 
with the kinetic energy operator H Q ; and 

Jn«(r,t) = Jdr{\i n (r\r),H'(r)}v a {Q a (r\t,o)}} (is) 

is a scattering operator in the Markovian approximation accounting for the effects 
of the collisions generated by the interactions with the surrounding media via an 
interaction Hamiltonian H', and where [cf. Eq. (A. 21)] 

V a {Q a (r | t, 0)} = [g a (r | t, 0)] Q /Tr{[g a (T \ t, 0)] a } , (19) 

that is, the escort probability that the unconventional statistics requires to be 
used in the calculation of averages values (cf. Eq. (18) in I). To solve the system 
of Eqs. (16) and (17), we need in Eq. (18) to express the right-hand side in 
terms of the basic variables. The scattering operator takes in general the form of 
the kind that is present in the relaxation-time approach, as for example shown in 
Refs. [23-25], namely 

J na (r,t) = -I n (r, t)/r Ia , (20) 

where Tj a is the momentum relaxation time (see Ref. [26], where it is presented 
an analysis of diffusion in the photoinjected plasma in semiconductors dealt with 
in the conventional statistical mechanics once the sufficiency condition is verified). 
Transforming Fourier in time Eq. (17) we have, after using Eq. (20), that 

(l + iUT Ia )I n (T,u) = -T Ia V-lW(T,u) , (21) 

which in the limit of small frequency, meaning UTj a <C 1, becomes, after back 
transforming to the time coordinate, 

I B (r,0 = -r /a V./a(r,0 , (22) 
and then, after using Eq. (22), Eq. (16) becomes 

£n(r,f) = r /a V./H(r,f) . (23) 



In order to close Eq. (23) we need to express the second-order flux in terms 
of the basic variables, n and I n , and after some calculus we find that (cf. Eq. 
(A. 20) in Appendix A where some additional comments and considerations are 
presented) 

V • (r, t) = £ a V rC- (r, t) , (24) 
where power index r y a has the a-dependent expression 

5 — 3a . . 

la = , (25) 
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where the values of a are restricted to the interval 1 < a < |, (see Appendix A) 
and consequently it follows the "anomalous" diffusion equation as given by Eq. 
(9), where D = C, a Tj a , once we take £ a and rj a as varying slowly in space and 
time, neglecting then their dependence on (r,t). 

Let us analyze this question of anomalous cyclic voltammetry in fractal elec- 
trodes in terms of the previous results. The solution of Eq. (9), for movement 
in one dimension, say, in re-direction normal to the electrode surface, is given, for 
< 7 < 1, by [27] 



n ( x , t) = b a r^ [a 2 + x 2 t~ 2 ^] ' 
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where a and b a are constants and /j a = {j a + 1)~ = | (3 — a) / (2 — a) . Hence, 
taking into account that the current, as noticed, is proportional to the arriving 
charges, and, once fi a > 0, admitting that x 2 t~ 2 ^ a <^ a 2 , there follows that 
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Finally, taking into account that the applied field is e (t) = e Q + v t (where v is the 
scanning velocity), and then t — (e — e a ) /v, Eq. (27) leads to the potential law 

/ ~ v* , (28) 

which stands for the peak value in the experiment and where £ is /i a , that is 
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(29) 



then it has been proved the empirical law, with the power £ expressed in terms of 
the infoentropic index a; the latter is then determined from the experimental value 
of £, which should be contained in the interval 0.5 < £ < 1, once 1 < a < 5/3. 



In Fig. 2 are shown the experimental results linking the values of the peak 
current for varying values of the scanning velocity v, from which can be derived 
the values of £ for each v, and then those of a through Eq. (29). Hence, keeping 
all other characteristics of the experiment fixed, the entropic index is dependent 
on the experimental protocol, in this case on v. From Fig. 2 and Eq. (25), it 
can be derived (similarly to what was done in the previous section) a "path to 
sufficiency", given in this case by 

a (v) ~ , (30) 

w v + v 2 

where, by best fitting, v 1 ~ 235 ± 20 and v 2 ^ 145 ± 15, in mV/s. As we can 
see, as v largely increases a goes to 1 - when one recovers the conventional result, 
that is, the standard Fick's law, indicating that the Criterion of Sufficiency in 
the conventional Boltzmann-Gibbs formalism is satisfactorily verified - while for 
small values of v, a tends to the lower bound value i>i/i>2 (~ 1.6 ± 0.15 in the 
experiment considered above). We can interpret this as a consequence that a 
large rate of change of the electric field leads to a very rapid transit of the charge 
through the roughened layer at the boundary giving no time for the movement 
tending to get adapted to the morphology of the thin film boundary. 

In the present type of experiments in electrochemistry we can derive the same 
conclusion as the one of the previous section in the experiments in semiconductor 
physics: the infoentropic index a is not an universal one for a given system but 
depends on the correct dynamics, geometry and size as well as the characteristics 
and influence of the boundary conditions, the macroscopic (thermodynamic) state 
of the system, and the experimental protocol. 

Two other interesting problems belong to this question of an "anomalous" 
diffusion-advection in polymer solutions and "anomalous" current in ionic-conducting 
glasses [28]. In the first case, macromolecules under flow, one faces the difficulty in 
the description resulting from the self-similarity in a type of average fractal struc- 
ture (what we have called "Jackson Pollock Effect" , in view of the anology with his 
painting with the dripping method, leading to fractal structures [29]). Without 
entering into details, given elsewhere [28], the "anomalous" diffusion-advection 
equation has the form 

tgn (r, t) - D a V 2 n^ (r, t) = -r a V ■ V • (n (r, t) [v (r,t) v (r,t)]) , (31) 
where r y a is the one of Eq. (25), D a of Eq. (9). 



In the second case, ionic-conducting glasses, the difficulty in dealing with 
them consists into the presence of long-range correlations in space, mediated by 
Coulomb interaction, having scaling characteristics. It leads to the so-called Curie- 
von Schwaitler current 

J{t)~r v , (32) 

that is, following a power law in time. Again, resorting to USM (in Renyi ap- 
proach) while ignoring the scaling characteristic, we arrive at such power law, with 
power v related to the infoentropic index by the expression 

-r^, ■ (33) 

(a - 1) 

4. IDEAL GAS IN INSUFFICIENT DESCRIPTION 

As additional illustrations we consider, first, an ideal gas of N particles in a 
volume V, but in the thermodynamic limit and we call ij and pj the position 
and velocity of the j-ih particle, with H (T) = ^p|/2m being the Hamiltonian. 

3 

The most general description of this system in any circumstances is in terms of 
all the observables of the system, that is from the knowledge of the single-particle 
dynamical function [5], which in classical mechanics is given by 

N 

n 1 (r,p|r) = ^5(r-r,) * (p - p,-) . (34) 
j'=i 

The statistical operator depends on it and, in NESEF (see I), on the associated 
Lagrange multiplier which we call ip 1 (r,p;t), determined by the informational 
constraints consisting in the average values 

m (r, p; i) = J dT ni (r, p | T) g(T\t,0) , (35) 

and we do have in the conventional formalism 

g(T\t,0) = Z^^exp!^- J d 3 r d 3 p Vl (r,p;t) ni(r,p|r)| , (36) 

for the auxiliary ("instantaneous frozen quasiequilibrium" ) statistical operator, 
with the nonequilibrium statistical operator given in terms of it by Eqs (1.14) to 
(1.17) in I. 



As shown elsewhere [5, 22] we can alternatively write the statistical operator 
in the form of a generalized grand-canonical one, the associated auxiliary operator 
being 

g GC (r I t,0) = Z~ GC (t) exp | - J d 3 r d 3 p [(3 (r, t) h (r) + A (r, t) n (r) + 

+i/ h (r, i) • I h (r) + v n (r, t) • I n (r) 

+ E K ( r > *) ® ^ ( r ) + ^ ( r > *) ® * (r)] ] } (37) 
In this Eq. (37) are present the quantities (generalized fluxes) 

/M ( r |r) = f d 3 p «W (p) nx (r, p | T) = ^ / rf 3 p «W (p) n„- (r, p | T,) (38) 

with r = 0, 1, 2, ... and 7 = h or n, and where 

n y (r, p I Tj) = 5(r- r 3 ) S (p - Pj ) (39) 

is the individual one-particle dynamical operator and Tj indicates the one-particle 
phase point (ij,Pj), and 

«M (p) = [u (p) ... (r - t^mes) ...u (p)] , (40) 

«M (P) = (p 2 /2m) U W (P) • (41) 

Moreover, in Eq. (40) [u (p) ... (r — times) ...u (p)] indicates the tensorial product 
of r-times the generating velocity 

u(p)=p/m , (42) 

where r = stands for the densities of energy, h (7 = h), and particles, n (7 = n), 
and their fluxes of all orders, i.e. the vectorial ones (currents) for r — 1, and 
the higher order (tensorial) ones for r > 2; moreover, we have designated the 
Lagrange multipliers by (3, A, u 1 ,F^ \ and we recall that ® stands for fully con- 
tracted product of tensors. For an alternative derivation of the auxiliary operator 
of Eq. (37), through the use of the auxiliary one of Eq. (36), see Appendix B. 



The average of the quantities of Eq. (38) over the nonequilibrium ensemble 
provides the macroscopic variables that are the basic ones for the construction of 
a nonlinear higher-order hydrodynamics (for example, in a linear approximation, 
see [24]). The single-particles are independent and then the statistical operator 
of Eq. (36) can be written as a product of distributions for each particle, namely 

N 

-Q(r\t,o) = l[-Q j (r j \t,o) , (43) 

with 

- Qj (T j | f , 0) = Zr l (t)exp ^- J d 3 r d 3 p Vl (r,p;i) h Xj (r,p | r,)| , (44) 

where h\j is the one-particle dynamical operator of Eq. (39), Zj ensures the nor- 
malization, tp 1 is the corresponding Lagrange multiplier. In the grand-canonical 
description we have that 

g 1GC (T j | t, 0) = Zr 1 (t) exp |- J d 3 r [/3 (r, t) Ah, (r) + A (r, t) An, (r) + 
+u h (r, t) • Ai ft i (r) + v n (r, i) • Ai nl (r) + 

+ E K 1 (r ' l) ® ^ (r) + F - (r ' t] ® AiS (r) ] 1 } ' (45) 



r>2 



where 



A^ 1 (r) = ^ 1 (r)-Tr{^ 1 (r)e 1 (r j |*,0)} , (46) 

/2 
^ 3 ^n(r,p|r,) , (47) 

and similarly for An, (r), etc. 

But let us now consider - as it was in the case of the preceding Section 3 - 
the situation in which we are forced to resort to a truncated description using, 
say, a finite small number s of fluxes (that is r = 0, 1,2, ...,s; s — 1 in Section 
3). Therefore Fisher's criterion of sufficiency is then violated if we proceed with 
the conventional statistics. Hence, if we want to analyze an experiment where 
these s fluxes are accessible, calculating expected values and response function 
depending on them, we should introduce USM. We proceed next to analyze such 



incomplete description of the hydrodynamics of the ideal gas in USM, resorting 
again to Renyi approach [2,3]. 

Consequently, for obtaining the one-particle statistical operator, g 1 (Tj | t, 0) 
one proceeds to maximize Renyi generating functional in terms of escort proba- 
bilities as discussed in I, that is [cf. Eq. (120)] 



nn 



(r,p;t) = J 



dTj ni(r,p | r,) V a {g^Tj \ t,0)} 



(48) 



which for a — 1, i.e. when the condition of sufficiency is satisfied, goes over the 
conventional expression, as given by Eq. (1.5). 

The heterotypical probability distribution is built in terms of the resulting 
auxiliary one, the latter, in the grand-canonical description, given in this case by 



Qi a (r, |t,o) 



1 + (a - 1) J d 3 r d 3 p [P a (r, t) Ahi (r) + A a (r, t) Arii (r) + 



Via (*) 

+v ha (r, t) ■ AI W (r) + v na (r, t) ■ AI nl (r) 



i 

l-a 



r=2 



(49) 



where A/ii, Ani, etc., are those defined after Eq. (47). The distribution of Eq. 
(49) can be written in an alternative form using the redefinitions of the Lagrange 
multipliers given by 



and 
where 



(3 a (r, *) = /?« (r,f) /G a (t) 
A a (r,t) = A a (r,t)/G a (t) 

"ha — v hal G a (t) , 
Vnoi = v nal G a (t) 

Fl[l(r,t) = FV(T,t)/G a (t) 



(50) 
(51) 
(52) 
(53) 

(54) 



G a (t) = 1 - (a - 1) J d 3 r d 3 p [f3 a (r, t) h (r, t) + A a (r, t) n (r, t) + 



+i^ Q (r, t) ■ l h (r, t) + v na (r, t) • I n (r, f ) + 



+ 



E[4 r i(r,t)®4 rl (r,t) + ^ r 2(r,t)®/M(r,t)]] , (55) 



r=2 



to obtain that 



l + (a-l) y rf 3 r rf 3 p £ Q (r,p;r, \ t) 



where 



with 



J d 3 r d 3 p £ a (r, p; T 3 | t) = tl ja ( Tj , Pj \ t) 



(56) 



(57) 



Qja ( r j) Pj 



= I a fo, + /3 a (Tj, *) + "na (Tj, t) • ^ + 



m 



2m m 



(5f 



^ + E E ^ 

r=2 7=n or /i 

after using the definitions of the densities and fluxes of order r = 1 to s in terms of 
the individual single-particle dynamical function of Eq. (34), and the integrations 
in r and p are performed.. 

Consequently, the auxiliary probability density of Eq. (43) becomes the prod- 
uct of factors involving each particle individually, given in the Renyi's approach 
we used by 



TV 



- Qa (r 1 1, o) ~ f[ [i + (« - 1) n,-« (r 1 1)}^ 

7a \ l ) j =1 



(59) 



In section 3 we have applied these results, for s = 1, to derive an "anomalous" 
diffusion equation (for fitting the experimental data) in voltammetry measure- 
ments in microbatteries with fractal electrodes, and in continuation we present 
an additional application consisting in a study of hydrodynamic properties of an 
ideal gas of photons in black-body radiation, extending these results to deal now 
with a quantum system. 



5. RADIATION IN INSUFFICIENT DESCRIPTION 



Let us consider an ideal gas of photons in the presence of an uniform flux of 
energy (generated, for example, having different temperatures at both ends of 



the container). We look for the determination of the energy, what is done, on 
the one side using the conventional approach in a description that includes the 
energy and the energy flux as basic variables and, on the other side, using a quite 
incomplete description which includes only energy, but dealt with in USM. We 

indicated by = c k the photon frequency-dispersion relation, and by (^Cj 

the annihilation (creation) operators in states |k). The energy and flux of energy 
dynamical operators are given by 

H = ^2 ala k 

k 

l h = ^ V k w k ala k . (61) 

k 

The auxiliary ("instantaneous frozen quasiequilibrium" ) statistical operator 
including both dynamical variables in the conventional formalism is 

g (t, 0) = Z' 1 (t) exp { -F h (t) H-v h (t)-I h } , (62) 

where F h and F h are the associated Lagrange multipliers (intensive nonequilibrium 
thermodynamical variable), with the nonequilibrium statistical operator given, in 
terms of this g, by Eq. (A. 14) to (A. 17). The averages of the basic dynamical 
variables can be calculated to obtain respectively that 

E(t)=^2 feuk ^k (t) , (63) 

k 

I h (t) = J2^ V ^k ^k (t) , (64) 

k 

where 

z/ k (t) = [exp {F h (t) fkJ k + v h (t) ■ ftu k V k ^k} - l] _i (65) 

Considering the presence of a weak flux, introducing an expansion in (t) 
keeping only the first nonnull contribution, we obtain the equations of state (i.e. 
relation between the basic nonequilibrium thermodynamic variables and the in- 
tensive one) given by 

i E(t)=aF- 4 (t)+a s (t) \u h (t)\ 2 , (66) 



(60) 



1 

V 



I h (t) = b (t) u h (t) 



(67) 



(68) 
(69) 
(70) 



where 

a = it 2 /15H 3 c 3 

a s (t) = 2tt 2 /9h 3 cF% {t) 

b (t) = -A7T 2 /A5h 3 cFl (t) 

as shown in Appendix C. It has been used the relation Wk = c k, and we notice 
that in equilibrium v h — and F h = [3 = [/c^T]" 1 , and then Eq. (66) becomes 
Stefan-Boltzmann law. 

Using Eq. (67) in Eq. (66) we can write for the energy that 

2 



v E(t) -m 



+ 



15 



V 



(71) 



Let us next consider the most insufficient description, done in terms of only 
the energy, but using USM in Renyi's approach, when the auxiliary statistical 
operator is given by 



Qa (*,0) = 



1 



Va (t) 



1 + («-!) AH 



(72) 



with AH = H — E (t), and we recall that in the calculation of average values is 
used the escort probability expressed in terms of g a [cf. Eq. (19)]. On the given 
assumption that the flux is weak, one should expect that a is close to 1, we write 
a = 1 + e (or e = a — 1) and we perform the calculations retaining only terms up 
to first order in e. As shown in Appendix C, after some calculation we find that 



where 



E(t) = E (t)-(a-l)<p h (t) a 2 E (t) 
E (t) = (t) = 

k 

= VT(4) C(4)/2tt 2 (Hc) 3 <pi(t) = 



(73) 



(74) 



with a given in Eq. (68), Y and ( being Gamma and Riemann functions respec- 
tively, which follow after using that 



A/k (*) = [exp {if h (t) hw k } - 1] 



(75) 



and a\ (t) is the fluctuation of energy 



4W = E(N 2 AT k (t) [l+A/- k (t)] 



AT k (t) 



4(a-l) [VT(4) C(4)/2vr 2 (fc) 3 y^)] 



and then 



ay?" 4 [l-4(a-l)] 



(76) 



(77) 



Equating Eqs. (71) and (77) there follows an equation for the infoentropic index 
a in terms of the variables which characterize the macroscopic state of the system, 
namely 



4 [l-4(a-l)] 



15 ^ 

8 ac 2 



I 

V 



1 2 



(78) 



But in the noted condition of weak flux we admit that i*X ~ f3 Q and ~ /3 
with, we recall, e = a — 1, and where /3 G is the reciprocal of the average temperature 
of the gas (its gradient being small). Hence, using these results in Eq. (78) we 
arrive at the expression for a given approximately by 



a 



1-C 



where 



15 J3l_ 
8 a 2 c 2 

,2 



(79) 



(80) 



It can be noticed that C 1 ~ (a/3 G 4 c) is something like the square of a flux 
of energy composed of the energy density of the radiation, a(3~ A , traversing at the 
speed of light, while we should expect / to be composed of something like the 
density of energy traversing at a speed determined by the gradient of tempera- 
ture, and then, in fact, we do have that for the ideal gas of photons (black-body 
radiation) works quite well meaning that the insufficiency of characterization can 
be ignored. 



6. IDEAL GAS IN A FINITE BOX 



We consider in this section other observation concerning the ideal gas, namely, that 
as a rule work in the thermodynamic limit is usually taken, what implies an infinite 



volume, but real systems are finite [30]. Then in the case of an ideal quantum 
gas in a finite box of volume V and area A, summation over the states k = 



n. 



(ir/L z )] (where n x ,n y ,n z are 1,2, etc. and L x ,L y ,L z 



are the size of the sides of the box) can be exactly replaced by an integral only 
in the thermodynamic limit. Otherwise we do have corrections depending on the 
area of the surface, and in the calculation of the partition function, internal energy, 
etc., there appears nonextensive terms. 

According to Pathria's textbook [31], and see also [32], the density of states of 
the ideal gas in the finite box is given by 



9(e) = V 



2vr (m/2n 2 h 2 f 2 e 1 ' 2 - ~ (m/nti 2 ) + 

8 V 



(81) 



where V is the volume and A the area, and dots stand for contributions quadratic 
and of higher powers in (A/V); once we consider a finite but large box it is kept 
in what follows only the first-order contribution in (A/V). 

Taking the equilibrium grand-canonical ensemble with temperature T and 
chemical potential /i, using the density of states of Eq. (81) and taking for sim- 
plicity the statistically nondegenerate condition, it follows that (see Appendix 
D) 



E 
N 



\k B T 



1 - - 



1 A\ T 



\T 3 exp{fi/k B T } 



for the density of energy, and 



N 
V 



1 - 



1 A\i 



4 V 



A T 3 exp {n/ k B T } 



(82) 



(83) 



for the density of particles; At is the mean thermal de Broglie wavelength, namely, 
= h 2 /rnksT. Hence, the energy per particle is given by 



1 + 



1 AX T 
12~V~ 



(84) 



Evidently, in the thermodynamic limit (infinite box) it follows the standard result 
E = (3/2) Nk B T (we recall that we took the nondegenerate condition). 

Next we admit to be in the condition of insufficiency consisting in ignoring the 
finite size of the box, and we proceed with the calculations in the thermodynamics 
limit but resorting - to patch the limitation thus introduced - to Unconventional 
Statistical Mechanics. 



Using USM in Renyi's approach, in the insufficient condition, namely, using 
the thermodynamic limit (infinite-size box), that is using in the calculations the 
escort probability [cf. Eq. (19)] in terms of the auxiliary Renyi's heterotypical 
distribution 



Q a = — {l + (a - 1) k (H - (if)) + F n ( 



N -IN 



i 

^ l-a 



(85) 



and under the expected condition of a near 1 (i.e. \a — 1| <C 1 for a very large 
box it follows that 

^-^-(a-^^al + F,^] , (86) 

N = N - (a - 1) [F n + F fe <7^] , (87) 

where a E and <7^ are the energy and particle-number correlation functions, and 
a 2 EN and a 2 NE the cross-correlation functions (see Appendix D, for details) and 
where E Q = (3N /2Fh) and iV are the first contribution in the expansion around 
a — 1. 

Finally, we can write (see Appendix D) 



E 3^ , 
~ -F~ 

N 2 h 



l + i(a-l) 



where, evidently, in the thermodynamic limit a approaches 1 as it should. More- 
over, from Eq. (88) we can clearly see that the infoentropic index a depends on 
the system dynamics, its geometry and size, and the thermodynamic state. 

Further applications of USM are available in the literature on the subject, and 
we can mention the interesting cases of its use on dealing with hydrodynamic 
turbulence and collider physics in [33] and [34] respectively. 



7. FINAL REMARKS 



We have illustrated in this paper the use of what we have called Unconventional 
Statistical Mechanics, an auxiliary formalism which can provide a theoretical ap- 
proach to situations when the conventional, and well established Boltzmann-Gibbs 
statistical mechanics has its use impaired because of a lack of a proper knowledge 
(for the problem in hands) of the characterization of the system and its dynamics 
on the part of the researcher (the Fisher's criterion of sufficiency is not satisfied 
as discussed in [1]). 



We have here dealt with "anomalous" luminescence from quantum wells in 
semiconductor heterostructures, where the failure of sufficiency resides at a mi- 
croscopic level, in that one does not know the proper quantum mechanical states 
in the thin (nanometric scale) quantum well. This is the result that we cannot 
precisely solve Schrodinger equation for the carriers because a failure to impose 
boundary conditions on the fractal-like morphological structure of the boundaries, 
to which we do not have easy access; evidently the wave function is largely affected 
by the roughness of the frontiers in the nanometric structure. As shown, as the 
width of the quantum well increases, the "anomaly" tends to disappear. This is 
an example in the realm of semiconductor physics. 

A second illustration, the one in section 3, is in the area of electro-physico- 
chemistry, namely the behavior of thin electrodes with a fractal-like morphology 
in microbatteries. An analysis of experiments of the so-called cyclic voltamme- 
try has been performed, whose results are a consequence of the occurrence of an 
"anomalous" diffusion of charges from the electrolyte. In this case the failure 
of sufficiency resides at a macroscopic, or, better to say, mesoscopic level, when 
classical Onsagerian hydrodynamics is used instead of the higher-order hydrody- 
namics that the problem requires. As described the "anomaly" depends on the 
experimental protocol, and tends to disappear as the rate of charge transfer is 
increased. 

The third illustration, see section 4, is a purely theoretical analysis concerning 
"anomalous" statistics of an ideal gas. As shown, once the ideal gas is perfectly 
and completely characterized in terms of the one-particle density function (Dirac- 
Landau-Wigner single-particle dynamical operator in the quantum case), the only 
possible description is the conventional one. But, if instead we work with trun- 
cated sets of linear combinations of the one-particle density, viz. the densities of 
particles and energy and the set of their fluxes up to a certain order, in a par- 
tial generalized grand-canonical description (and then we do have insufficiency 
of description), to overcome the fact that we are not complying with the crite- 
rion of sufficiency (i.e. to calculate properties) unconventional statistics can be 
introduced, at the price of having undetermined parameter (s). In terms of them 
one obtains "anomalous" hydrodynamic laws, for example, in the lowest order 
approximation, "anomalous" diffusion, with a power law, in the subsequent order 
an "anomalous" Maxwell-Cattaneo equation (damped waves), and so on. As the 
order used in the higher-order hydrodynamics increases the "anomaly" tends to 
disappear because one is approaching a description equivalent to the one in terms 
of the single-particle density matrix. 



It can be noticed that we have always written the word anomalous within 
quotation marks: This is so because there is nothing anomalous in the physical 
laws governing the system: what is "anomalous" is the result obtained using 
the unconventional approach when, as noticed, we are unable to comply with 
the criterion of sufficiency in the use of the proper Boltzmann-Gibbs statistics. 
Moreover, concerning the, what we have called, "path to sufficiency" evidenced 
in the three cases we presented, it is tempting to conjecture that it may be a 
general characteristic of unconventional statistical mechanics, in any of the cases 
when different infoentropic-index-dependent structured informational entropies 
are used. 

Finally, we add a couple of comments concerning other possible situations that 
can be analyzed in terms of USM, and its connection with the criteria of efficiency 
and sufficiency 

Application of a USM resorting to Havrda-Charvat structural infoentropy to 
the question of the evolution of a nonequilibrium temperature-like variable, say 
T* (t), in the case of a system with long-range spatial correlations [35] shows that 
in the approach to equilibrium is present a long and near-stationary plateaux 
above the reservoir temperature which is to be attained in the long range once 
final thermal equilibrium follows. The calculations are not related (compared) to 
experimental results, and so we do not know if they have any meaning, except to 
leave a suggestion of theoretical results in search of an experiment. If the exper- 
iment shall show in fact that there are differences in relation to the conventional 
calculation ("anomalous" results) one should look where Fisher's criteria of effi- 
ciency and sufficiency is not satisfied (if it is the case, as already noticed should 
be related to improper handling of the long-range correlations; see below the case 
of the one-electron transistor). There exists a situation in the physics of semicon- 
ductors where is also present a near stationary state in the way to equilibrium 
associated to a slowing-down of the relaxation processes. In the case of the out 
of equilibrium photoinjected double plasma in semiconductors (see for example 
Chapter 6 in [5]) what is called the carriers' quasitemperature T* (t) follows a 
path to equilibrium which strongly depends on the experimental protocol [36]. 
Thermal equilibrium with the reservoir follows very rapidly (tens to a hundred 
picoseconds) when the exciting laser pulse is short (pico- to subpico-second scale), 
but is slow and presenting a long near plateaux when the exciting lase pulse is long 
(tens of picoseconds duration). Such behavior in the latter case is not properly 
described by theory in the conventional formalism but could, we think, as in [35] 
be obtained resorting to USM with an appropriate adjustment of the infoentropic 



index involved in, for example, Havrda-Charvat or Renyi approaches. The original 
failure in this case consisted that in the treatment that was applied one does not 
comply with the criterion of sufficiency by not considering the facts responsible 
for the phenomenon to occur, namely, "phonon bottleneck" and ambipolar diffu- 
sion (the latter keeps decreasing the density of particles), which once incorporated 
in the description produces an excellent agreement of (conventional) theory and 
experiment, and we have sufficiency at work [36]. 

Another interesting case to be considered is the one of single-electron transis- 
tors [37], when long-range and strong (because not screened) Coulomb interaction 
leads to difficulties when dealt with in a simple way in the conventional approach. 
Again, the point is a failure of sufficiency: in this case is the proper charac- 
terization of the states of the system as a result of long-range interconnected 
correlations. When renormalization of the carrier states is introduced, sufficiency 
is restored and there follows an excellent agreement between (conventional) the- 
ory and experimental data [37]. We conjecture that instead of application of the 
renormalization group for, as said, restore sufficiency (i.e. a good mimical de- 
scription of the carriers' states), one could use the unrenormalized description in 
terms of USM (a point to be considered in the future). 

In the last two cases above it is also conjecturable that there may be present 
a kind of "path to sufficiency", of the like of those presented in the three other 
former examples. 

In a concluding remark, we can say, in summary, that the illustrations here 
presented of the application of the theory of the preceding article, allows for gain- 
ing a better perspective of USM, which appears as a useful and practical formalism 
for the macroscopic description of systems when the research cannot comply with 
criteria of efficiency and/ or sufficiency in the conventional, well established, and 
physically and logically sound Boltzmann-Gibbs theory. Basically it is a sophisti- 
cated fitting formalism to be used in the referred conditions, namely, when one 
does not have a proper access to the characteristics of the system that are relevant 
to determine the property of the system we are studying. 

It is worth noticing that in the illustrations we presented there appears a quite 
interesting connection between theory, experiment, and front-line technology. 
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Appendix A. "Anomalous" Diffusion 

Let us consider first the conventional case of diffusion, when the criterion of suf- 
ficiency is satisfied, meaning the several stringent restrictions its validity requires 
are met, namely, when are satisfied the conditions of local equilibrium, linear On- 
sager relations and symmetry laws, condition such that the motion is dominated 
by long wavelengths and very low frequencies contributions, and weak fluctua- 
tions. A specific criterion for validity is given in Ref . [38] , where it is considered 
a Brownian-like system composed of two ideal fluids in interaction between them. 
The continuity equation for the flux I n (Eq. (17) in [38]) after transforming Fourier 
in time, takes the form 

(l + iur nl ) I n (r,u)+r nl W ■ /| ] (r, cj) = , (A.l) 

that is Eq. (10) for r — 1, and using Eq. (20) for the collision integral and 
then r n i is the momentum relaxation time. But, at sufficiently low frequencies, 
meaning that wr„i can be neglected, we have the appropriate expression 

I n (r,u;)~-T nl V- /J?(r,cu) , (A.2) 

and a direct calculation tells us that 

V • /| ] (r, w) ~ (k B T/m) Vn (r, u) . (A.3) 

Replacing Eq. (A.3) in Eq. (A.2), and the latter in the conservation equation 
for the density (Eq. (16) in [38]), we obtain the usual Fick's diffusion equation 

^-n(r,t)+D V 2 n(r,f) = , (A.4) 

where D = \v\ h r nl , with |m v\ h = ^ksT; v t h is the thermal velocity and D the 
diffusion coefficient, with dimension cm 2 / sec. 

Let us now go over the unconventional treatment, which is required once one 
is looking forward for an adjustment of data on the basis of a description in 
terms of a diffusive movement, when this is not possible, as a consequence that 
diffusion in the microroughneessed region is governed by not too long wavelengths 



(up to the nanometric ones, i.e. 10~ 7 cm, while the limitation of the diffusive 
domain [38] is the order of D/v t h, say, typically 10~ 2 to 10~ 4 cm). A higher- 
order hydrodynamics [24] needs be introduced, but if the lower order description 
including only the density and its flux is kept, then we are not complying with 
the criterion of sufficiency and we need to introduce Unconventional Statistical 
Mechanics. Let us consider the auxiliary statistical operator which for this system 
of free particles can be approximated by a product of the statistical operator of 
the individual particles, namely (cf. Section 4) 



A? 



g(r\t,o) = ]lQ la (r 1 \t,o) 

3=1 



(A.5) 



where F 1 stands for the one-particle phase point (rj,pj) and in Renyi statistics 
we do have that 



Qia (r x |t,o) = 



z a (t) 



1 + («-!) / d 3 r / d 3 p (p la (r,p,t) An n (r,p,t) 



where 
with 



and 



An n (r, p, *) = n u (r, p | I\) - (nn (r, p | T 1 )) a 

nu (r, p | r\) = 5 (r - Tj) 5 (p - pj) 
(nn (r, p | rO) Q = JdT fin (r, p | T 1 ) Via {Qi a {Vi | t, 0)} 



P P 

ip la (r, p, t) = F na (r, t) + F ha (r, t) h F na (r, t) 



(A.6) 
(A.7) 

(A.8) 
(A.9) 



2m m 

are modified forms of the associated Lagrange multiplier which appear in Eq. 
(1.16) in I, where in the latter the change in space of the energy density has been 
desconsidered (i.e. we took Fh a (r,t) = f3), z ensures its normalization, and 

Via {Qia 1 1,0)} = [g la (r x 1 1, 0)] Q / J dTx [g la (ri I t, 0)] a , (A.10) 



is the accompanying escort probability (see I; [2]). 



Introducing the modified Lagrange multiplier 



Via (r,p,t) = Via (r,P,0/ 
we find that 



I - (a - 1 ) y d s r J d 3 p<p la (r, p, t) (n u (r, p | r\)) Q 

(A.H) 



ft a (r 1 1, o) = 



i 



z(0 



1 + (a - 1) / d 3 r / rf 3 p £> la (r, p, i) n u (r, p | r x ) 



(A.12) 

with Z (t) ensuring the normalization condition. 

Using Eq. (A. 5) to Eq. (A.12), after some lengthy but straightforward calcu- 
lations, it follows for the energy density that 



h (r, t) = u (r, t) + n (r, t) -mv 2 a (r, t) 



(A.13) 



i.e., composed of the energy associated to the drift movement (the last term) and 
the internal energy density 



where 



C a (r,t) 



u(r,t) 



2ixN 



3 C a (r,t) 
5 - 3a (3 a (r, t) 



n 1a (r, t) 



(A.14) 



(a-l)'Z(t) 



2m 



-B 



a 



P a (r,t) \2'a-l 2 



2("-l) 
a-3 



(A.15) 



B(v,x) is the Beta function, we have written F ha (r,t) = f3 a (r,t); F na (r,t) = 
mf3 a (r, t) v a (r, t) introducing a "drift velocity" field v a (r, t); moreover 



3a — 5 

lex 77 

a — 3 

and the values of a are restricted to the interval 



!<«<- 

Finally, the second order flux is given by 



(A.16) 



(A.17) 



I [ n ] M) = / 



d 3 p 



P P 

m mi 



(nu (r,p | r0) Q 



= n (r, t) [v (r, t) v (r, *)] + (r, t) 1® , (A.18) 

where l' 2 ' is the unit second order tensor, [...] stands for the tensorial product of 
vectors, and it can be noticed that 

P [2] (r, t)=m I [2] (r, t) - m n (r, t) [v (r, t) v (r, *)] (A. 19) 

is the pressure tensor of classical hydrodynamics. Neglecting the terms quadratic 
in the drift velocity, combining the above equations we obtain that 

4 2l M)=£„ a nMr,t) , (A.20) 

where 

£ 2 5 C °( r '*) (A21) 

W 3m5-3a p a ( r ,t) ' 1 ' ) 

is the quantity present in Eq. (24). 

Appendix B. The Grand- Canonical Probability Distribution 

We introduce now the Taylor series expansion in p in the Lagrange multiplier of 
Eq. (36), namely 

<p x (r, p; f) = T (r, t) + F (r, t) ■ p + F [2] (r, t) ® [pp] + ... , (B.l) 

where T — (p 1 with p = 0, [F]. = d(p 1 /dp i (i = x,y,z) with p = 0, is the 
tensor of components Tij = d 2 (f 1 /dp i dpj with p = 0, and [pp] is the tensorial 
product of the vectors p and ® stands for fully contracted tensorial product, etc. 
Hence, for the quantity Zqc (t) which ensures the normalization of the distribution 
we do have that 

2(r | t) = ^-^y exp |- J d 3 r d 3 p \T a (r, t) + F a (r, t) • p+ 

+^(r,t)®[pp] + ...] nx^plT)} . (B.2) 
Rewriting the Lagrange multipliers T as 

F(T,t) = A(r,t) , (B.3) 

F (r, t) ■ p = v n (r, • - + u h (r, t) ■ ^ , (B.4) 

m 2m m 



(r, t) ® [p... (r - times) ...p] = (r, t) ® «M (p) + if ] (r, f) ® «W ( p ) , 

(B.5) 

we recover Eq. (37) [we recall that the Lagrange multipliers for energy and density 
are not independent but related by the law as shown in Ref. [39]. 



Appendix C. Radiation Under Flow 

Using the auxiliary distribution of Eq. (62) - which, we recall produces for the 
basic variables and only the basic variables, in this case H and 1^, the same 
average value as the proper nonequilibrium distribution, cf. preceding article -, 
it follows that 

E (t) = Tr jx> k a[a k g (t, «)} = E ^ k Tr {a[a k g (t, 0)} , (C.l) 

N(t)=Trl^2 a W g(t,0)J=X>{ a[a k g(t,0)} , (C.2) 
and after calculation 

Tr{4a k S(*,0)} = i/ k (t) (C.3) 

with z/ k (t) given by Eq. (65) [a shifted-Planck population]. Considering the case 
of a weak flux, we can write 

u k (t) = [exp {F h (t) fko k } - I}' 1 + ... (C.4) 

where ... stands for an expansion of Eq. (65) up to second order in v^. Using Eq. 
(C.4) in Eq. (C.l) and (C.2) one obtains the expression of Eqs. (66) and (67). 

Let us now consider the insufficient condition in which the flux is ignored, and 
we resort to USM in Renyi's approach, introducing the heterotypical distribution 
of Eq. (72), is the one of Eq. (73) where 

E o {t) = (H) = Tr[H Qa {t^)} , (C.5) 

a 2 e (t) = (HH) o -(h) 2 o = Tt{(h- (h)) 2 ~g a (t, 0) j , (C.6) 

after performing an expansion around a — 1 in the exponent of the escort proba- 
bility (cf. Appendix A in I) and retaining the first contribution assuming a near 



1; it can be noticed that these average values are dependent on the infoentropic 
index a. Moreover, 

E Q (t) = Tr I fkj k ala k g (t, 0) > = ^ hj^Tr ja k a k £? Q (t, 0)| = ^ hu^Af^ (t) , 

Ik J k k 

(C.7) 

where jV k (t) is given in Eq. (75). 

Appendix D. Ideal Gas in a Finite Box 

In the conventional approach using the grand-canonical ensemble we do have the 
well known results that the average energy is given by 

oo 

E = Wk„ = Jdeg(e) ef (e) , (D.l) 

k<X 

and the average number of particles by 

oo 

iV = ^/ k „ = J deg(e) /(e) , (D.2) 

k<T 

where 

/( e ) = [exp{/?( e -/i)} + l]- 1 , (D.3) 

with (3 = 1/ksT and \x is the chemical potential, and = H 2 k 2 /2m. 

Using for density of states g (e) the expression of Eq. (81) in Eqs. (D.l) and 
(D.2), after taking the statistically nondegenerate limit, i.e. neglecting in Eq. 
(D.3) 1 in comparison with the exponential, there follow the expressions of Eqs. 
(82) and (83) in the text. 

We consider now the unconventional approach, introducing the insufficiency 
in Fisher's sense of taking for g (e) only the first term on the right of Eq. (81) - 
meaning the expression corresponding to the thermodynamic limit -, and resorting 
to Renyi's approach, that is, using the statistical operator of Eq. (85) and in the 
calculation of averages the escort probability in term of it, it follows that the 
average value of energy is given by 



Ea(H)-{a-l)F h 



BH \-{")] 



[a 



1 > F -[< M >. 

(D.4) 



and for the average value of the particle number 



N~(N)-(a-l)F n 
where 



-{a-l)F h 



(...)= Tr{...g a } 



(™) 

(D.5) 
(D.6) 

and it has been taken a is near 1, i.e. the case of a finite but long box. 

Performing the calculations, a lengthy but straightforward task, it follows that 

E = \KF^ -{a -I) (jN F h ^ + \N F n F-^ , (D.7) 



N 



N -(a-l)N (- + F n 



(D.8) 



It should be noticed that the averages of the kind of Eq. (D.5) present in these 
expression are dependent on a. If, in a first approximation we approximate these 
average by taking a ~ 1, from Eq. (D.7) and (D.8) there follows Eq. (87) and 
consequently the expression for the infoentropic index a as given by Eq. (88). 
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Values of the infoentropic index a and kinetic temperature fitting the 
luminescence spectrum as shown in Figure 1. 
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Figure 1: Comparison of theory (full line) and experimental results (filled geo- 
metric figures) in the luminescence spectra (high frequency side or Shockley- 
Roosbroecke domain purely dependent on the carriers' dynamics), for several 
samples under identical processes of growth, but with different length of the 
quantum well. The corresponding values of the infoentropic parameter a 
and the kinetic temperature are given in Table I. 
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gure 2: Logarithmic plot of the current peak value with the scanning velocity 
v of the applied electric field, in experiments of cyclic voltammetry in mi- 
crobatteries with nanometric thin film fractal-like electrodes. The tangent 
at each point gives the value of the index £ of Eq. (29); the lower and upper 
straight lines are approximate fittings in regions where £ varies smoothly [cf. 
Eq. (31)]. 



